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Electron injection process at high Mach number collisionless quasi- 
perpendicular shock waves is investigated by means of one- dimensional elec- 



tromagnetic particle-in-cell simulations. We find that energetic electrons are 
generated through the following two steps: (1) electrons are accelerated nearly 
perpendicular to the local magnetic field by shock surfing acceleration at the 
leading edge of the shock transition region. (2) the preaccelerated electrons are 



further accelerated by shock drift acceleration. As a result, energetic electrons 
are preferentially reflected back to the upstream. Shock surfing acceleration pro- 



Oh" vides sufficient energy required for the reflection. Therefore, it is important not 

only for the energization process by itself, but also for triggering the secondary 
acceleration process. We also present a theoretical model of the two-step acceler- 
ation mechanism based on the simulation results, which can predict the injection 
efficiency for subsequent diffusive shock acceleration process. We show that the 
^ ■ injection efficiency obtained by the present model agrees well with the value ob- 

tained by Chandra X-ray observations of SN 1006. At typical supernova remnant 
shocks, energetic electrons injected by the present mechanism can self-generate 
upstream Alfven waves, which scatter the energetic electrons themselves. 

Subject headings: acceleration of particles — cosmic rays — plasmas — shock 
waves 



1. INTRODUCTION 

The origin of nonthermal emission observed from a variety of astrophysical objects is 
still a major unresolved issue of plasma astrophysics. These include supernova shocks, extra- 
galactic radio sources and active galactic nuclei. Among them, shocks of supernova remnants 



-2- 



(SNRs) are believed to be the most probable acceleration site of Galactic cosmic rays. Much 
theoretical and observational work has been devoted to particle acceleration processes around 
shock wave s. One of t he most widely applied theory is diffusive shock acceleration (DSA) 



theory (e.g. lBelllll978aU bl; iBlandford Ostrikerlll978l ). It has been very successful in provid- 



ing a natural explanation for the power law distributions of high energy particles observed 
in many astrophysical sources. This process utilizes Alfven wave turbulence as the particle 
scatterers. Under the assumption of elastic collision between particles and waves, energetic 
particles scattering back and forth across the shock front gain a net momentum because of 
the converging velocity fields. The theory was extended to include the finite shock size and 
the self-consistent wave excitation and applied to the Earth' s bow shock in order to account 
for the diffuse ion component observed in the upstream (e.g. lEichlerlll98ll ; lLeelll982l ). While 



in situ observations of energetic ions associated with int erplanetary shocks and planetary 



bow shocks can be well explained by the DSA theory (e.g. iScholer et al.lll980l ; iGosling et al. 
19811). energetic elect rons thought to be accelerated by DSA process are rarely observed 
( jshimada et al. 1999 ). On the other hand, there is no doubt about the existence of ultra- 
relativistic electrons which may be accelerated by DSA process at SNR shocks. We still have 
no clear consensus on what determines the electron acceleration efficiency The physics of 
electron acceleration at collisionless shocks is poorly understood so far. 

The well-known difficulty is that thermal electrons cannot be easily scattered by Alfven 
waves because of their small gyroradii. Injecti on from therma l pool to mildly relativistic 
energy by some other mechanisms is required. iLevinsonl (119921 ) has examined electron in- 
jection process at strictly parallel shocks. He considered the self-consistent excitation of the 
whistler waves with cosmic ray electrons by applying the standard quasi-linear theory. It 
was shown that the injection of low energy electrons by the self-generated whistlers may 
be possible when the Alfven Mach number exceeds 43/a//^, whe re (3 P is the ra t io of ther- 
mal electron pressure to magnetic pressure. On the other hand, iPapadopoulosi (119881 ) has 
taken a different approach. He considered an electron energization process within the shock 
transition region by strong plasma microinstabilities. It is well known that the reflection 
of upstream ions plays a dominant role in the structure of quasi-perpendicular shocks (e.g. 
Leroy et al.lll982l ). The reflected ions streaming relative t o the upstream plasm a could excite 
various plasma instabilities in the so-called foot region. IPapadopoulosi (119881 ) argued that 
the relative drift velocity exceeds the electron thermal vel ocity at high Ma ch number shocks 
and leads to the excitation of the Buneman instability (iBunemanl Il958l ). The Buneman 
instability gives rise to very rapid electron heating. As a result, the interaction between in- 
coming/ refl ected ions and the prehe a ted el ectrons permits the excitation of the ion acoustic 
instability. ICargill fc Papadopoulod (119881 ) demonstrated the idea of the electron energiza- 
tion process by using a hybrid simulation code where ions are treated as particles whereas 
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electrons are assumed to be a massless charge-neutralizing fluid. IShimada fc Hoshinol (120001 ) 
extended their studies to include the electron dynamics in the shock structure by using a 
particle-in-cell (PIC) code where both ions and electrons are treated as particles. They 
found that localized large amplitude electrostatic solitary waves (ESWs) are produced in 
the nonlinear stage of the Buneman instability. They also argued that the rapid electron 
heating and acceleration are involved with ESWs. Now, the acceleration mechanism associ- 
ated with ESWs is considered as electron shock surfing/surfatron acceleration (SSA) process 
( jKatsouleas fc Dawsonlll983[ h 



S SA mechanism for ions has been extensively studied by many authors (e.g. ISagdeev 



19661 ; iLee et al.lll996t IZank et al.lll996l ). The process utilizes the electrostatic shock poten- 
tial which is caused by inertia difference between ions and electrons in the shock transition 
region. Ions having energy smaller than the shock potential are reflected by the shock front 
and begin to gyrate around the upstream magnetic field. During their gyromotion in the 
upstream, they are accelerated parallel to the shock surface by the motional electric field. If 
the spatial scale of t he shock potentia l is small compared to the ion inertial length, multiple 



reflection can occur (jZank et al.lll996l ). In contrast to this, electrons cannot be reflected by 



the shock potential. The electric field directed to the downstream is required. iHoshinol (120011 ) 
argued that ESWs can play the similar role to the shock potential in ion shock surfing ac- 
celeration, because ESWs are associated with phase space electron holes (positively charged 
structures). Electrons trapped in ESWs can be accelerated by the motional electric field. 
This acceleration mechanism is very efficient, so that mildly relativistic electrons are gener- 



a considerable attention and investigated in detail by many authors (e.g. 


Dieckmann et al. 


2000a: 


McClements et al. 


2001; 


Hoshino & Shimada 


2002; 


Schmitz et al. 


2002a 


o), because 



it may provide a clue to the electron injection problem. 

In situ observations of the Earth's bow shock also evidence the existence of ESWs 
in the shock tra n sition region and associated electron heating (e.g 



Bale et al. 1998, 2002 



Hull et al.ll2006l ). lOkal (120051 ) has recently carried out a detailed investigation of the well- 
resolved high Mach number Earth's bow shock (Ma — 14) observed by the Geotail satellite. 
He showed a clear coincidence between nonthermal electrons and the appearance of the broad 
band electrostatic noise (BEN) which is considered to be a signature of ESWs, although a 
plausible acceleration mechanism remains unanswered. 

Most of previous studies concerning electron energization processes by microinstabilities 
at high Mach number shocks are restricted to the case of strictly perpendicular shock. On 
the other hand, a number of PIC simulati ons of quasi-perpendicular supercritical shock s 



have been conducted for many decades (e.g. iForslund et al.lll984t iLembege k, Savoinilll992l ). 
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However, the main applications of these studies were planetary bow shocks and interplan- 
etary shocks in the heliosphere. Here we extend these studies to very high Mach number 
quasi-perpendicular shocks. We find that SSA produces suprathermal electrons in the tran- 
sition region of quasi-perpendicular shocks. The electron energization via SSA occurs within 
a relatively narrow region at the leading edge of the shock transition region where strong 
electrostatic waves are observed as in the case of strictly perpendicular shocks. The differ- 
ence is that the preaccelerated electrons are further accel erated by the so-cal l ed fas t Fermi 
acceleration process which is proposed by IWul (119841 ) and lLeroy fc Mangeneyi (11984J ) . They 
considered a particle motion in the de Hoffman- Teller frame (HTF) where the motional elec- 
tric field vanishes. In that frame, a particle having sufficiently large energy is reflected by 
the shock which acts as a fast moving magnetic mirror. The mirror reflection is adiabatic 
process, provided that a particle gyroradius is much smaller than the shock thickness. There- 
fore, the energy of reflected particle is conserved in the HTF. Then, after the interaction with 
the shock, the particle momentum parallel to the magnetic field measured in the upstream 
frame is increased by Ap = 2mV\/ cos#b„, where m, V± and 9bu are respectively, particle 
mass, the upstream plasma bulk velocity and the shock angle. Since the momentum increase 
is proportional to the reciprocal of cos#e n , the acceleration becomes extremely efficient at 
nearly perpendicular shocks. While the acceleration efficiency increases with increasing the 
shock angle, the number of reflected particle rapidly decreases, because the initial energy 
required for the reflection increases. 



Later, iKrauss-Varban &: Wul (119891 ) has shown that the fast Fermi process in the HTF 
is equivalent to shock drift acceleration (SDA) in the normal incidence frame (NIF) where 
the upstream velocity is parallel to the shock normal. In the NIF, a particle gains its 
energy by drifting pa r allel to the upstream motional electric field direction. Note that 
Krauss-Varban fc Wul (119891 ) called the acceleration mechanism as gradient drift acceler- 
ation and distinguished it from SDA. Since SDA was usually based on the approximation 
that a particle gyroradius is large compared to the shock width and a particle has multiple 
interaction with the shock front. In this paper, however, we will not discriminate the differ- 
ence and simply call the acceleration mechanism as SDA, because the physical mechanism 
is the same (gradient drift provides the energy gain). 

The process has been extensively studied in order to account for observed energetic 
electr ons in the upstream of the Earth's bow s hocks (e . g. IKrauss-Varban et al.lll989l ; IVandas 
19951 ). However, detailed parametric survey of lVandas! (120011 ) has demonstrated quantitative 
discrepancies between the theoretical expectation and observations. He argued that the pro- 
cess should be modified by some other nonadiabatic processes such as pitch angle scattering 
in order to explain observations. In addition to this, the required energy for the reflection of 
thermal electrons at very high Mach number shocks is unrealistically high. This process by 
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itself cannot account for the observed nonthermal electrons at high Mach number shocks. 



Previous studies did not consider the effect of microturbulence in the shock transition 
region, which we show plays an important role in the generation of energetic electrons. In 
fact, SSA is a highly efficient acceleration mechanism, so that the preaccelerated electrons 
gain sufficient energy required for the reflection process. As a result, the preaccelerated 
electrons are further accelerated via SDA. This two-step acceleration mechanism is important 
for providing a seed population for subsequent DSA process. In fact, we will show that the 
reflected electron energy is large enough to be accelerated by DSA process when the Mach 
number is typical of SNR shocks. 

We propose a theoretical model of the two-step acceleration mechanism based on the 
simulation results. The present model well explains the observed injection efficiency and the 
energy density of cosmic ray el ectrons which were obtained by detailed analysis of Chandra 
X-ray observation of SN 1006 (IBamba et al.l 120031 ). Moreover, the present model predicts 
the shock angle dependence of the injection efficiency. This dependence again agrees well 
with th e shock angle constraint, which is required to account for the observation by the DSA 
theory (IBamba et al.ll2003l ; lYamazaki et al.ll2004j ). 



2. SIMULATION 

We study the dynamics of ions and electrons in the self-consistent shock structure by 
utilizing a one- dimensional electromagnetic PIC code where both ions and electrons are 
treated as particles. A high-speed plasma consisting of electrons and ions is injected from the 
left-hand boundary of a one-dimensional simulation system and travels toward the positive 
x. The plasma carries a uniform magnetic field B x and B z . At the right-hand boundary, 
the particles are specularly reflected. Then, a shock wave is formed and propagates in the 
negative x direction. The downstream bulk speed becomes zero on average in the simulation 
frame. Initially, there are 100 particles for each species in each computational cell. The 
grid size is comparable to Debye length and the simulation box consists of 51200 grids. The 
plasma parameters are as follows: The upstream plasma (3 e = /3j = 0.08 (/3j = SnnTj/B 2 ), 
where n, Tj, B are the density, temperature, and magnetic field strength, respectively. The 
ratio of the plasma frequency to the electron cyclotron frequency is u pe /fl ce = 20. In order 
to reduce the computational costs, the ratio of ion to electron mass mj/m e = 100 is used. 
The upstream Alfven speed becomes 5 x 10~ 3 c, where c is the light speed. We use the 
plasma injection four velocity of Uq = 5 x 10~ 2 c. The Alfven Mach number of the resultant 
shock wave is Ma — 15 in the shock rest frame. Several runs with different shock angles are 
conducted with keeping the upstream magnetic field strength unchanged. In this section, we 
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mainly discuss the results of a run with Qbu = 80°. 

Figure [T] shows an overall structure of the shock transition region at uj pe t = 12000 
(corresponding to VL ci t = 5.5). From the top panel, ion phase space diagram in (X, Ui X ), 
electron phase phase space diagram in (X, U ez ), (X, U ex ) and (X, K e ), the magnetic field B z 
and the electric field E x , respectively. The plasma four velocity is normalized to the injection 
velocity Uq. The electron kinetic energy K e = (p{ e — l)m e c 2 is normalized to the injection 
energy K e0 = (70 — l)m e c 2 , where 7 e and 70 are the Lorentz factors of each particle and the 
injection velocity, respectively. The magnetic field and the electric field are normalized to the 
z component of upstream magnetic field Bq z = Bo sin 6 Bn and the corresponding motional 
electric field E 0y = UoB 0z /j c, respectively. The spatial scale is given in units of the electron 
inertial length c/u pe in the upstream. The color of the phase space diagrams represents 
the logarithm of particle count in each bin. Note that the vertical scale of the fourth panel 
(electron energy spectra) is also shown in a logarithmic scale. 

The basic structure of the shock transition region is simila r to those obtained by pre 



vious simulation studies o f strictly perpendicu lar shocks (e.g. IShimada fc Hoshinol 12000 



Hoshino fc Shimadal I2OO2I : ISchmitz et all l2002al lbh. We can find two distinct ion compo 



nents, the incoming and the reflected ions. Strong electrostatic turbulence is observed at 
the leading edge of the shock transition region (280 < X < 290), where the relative drift 
velocity between the incoming electrons and the reflected ions becomes maximum. These 
electrostatic waves are excited by the Buneman instability. Strong energization of the up- 
stream electrons coincides with the turbulent electrostatic waves as in the case of strictly 
perpendicular shocks. The heating and acceleration of electrons due to the turbulence occur 
on a very fast time scale. If we look at the deeper inside the shock transition region, the pre- 
heated electrons trigger the ion acoustic instability and the associated heating of incoming 
ions is evident in the top panel. The maximum growth rate of the ion acoustic instability is 
> lOQ^ 1 , which we obtain by using the standard linear dispersion analysis. The growth rate 
is large enough for the instability to develop within the shock transition region. Electrons 
are slowly heated up by the ion acoustic waves and the adiabatic heating process. 

In addition to these features, which are common to strictly perpendicular shocks, we 
can clearly find energetic electrons streaming away from the shock front along the magnetic 
field. These parallel escaping energetic electrons can be seen in the second panel, which 
represents the z component of electron four velocity. Note that the z component of velocity 
is almost parallel to the magnetic field. Figure [2] displays the distribution functions of 
electrons in (U X ,U Z ) plane taken at four different locations around the shock transition 
region. Just before the shock front (a), energetic electrons are observed as a distinct beam 
component. The typical beam drift velocity parallel to the magnetic field is u\\/Uq ~ 15. 
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Fig. 1. — Overall shock structure of high Mach number quasi-perpendicular shock 
80°). Color represents the logarithm of particle count in each bin. 
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Fig. 2. — Distribution functions of electron in (U X ,U Z ) plane. Each panel shows the dis- 
tribution function taken at (a) 260 < X < 280, (b) 280 < X < 300, (c) 300 < X < 320, 
(d) 320 < X < 340, respectively. The solid line represents the direction of the averaged 
magnetic field. Color represents the logarithm of particle count in each bin. 
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Fig. 3. — Energy spectra of electrons obtained around the shock transition region. Electron 
energy is normalized to the upstream bulk energy. 
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Since the electron beam excites Langmuir waves in the upstream, the cold upstream electron 
component is slightly modified from the far upstream condition. With increasing penetration 
into the shock, the incoming electrons are accelerated/heated up mainly perpendicular to 
the magnetic field, while the parallel drift velocity of the energetic electron beam decreases. 
Eventually, these two components merge into a single, but non-Maxwellian distribution (d). 
Figure [3] shows energy spectra of electrons. In the transition region (280 < X < 340), 
the middle energy range (10 < K e /K e0 < 100) of the spectra can be approximated by the 
power law with indices of 3 — 4. On the other hand, humps are found in the high energy 
part (K e /K e Q > 100) of the spectra. We observe the humps of the energy spectra only in 
the upstream or at the leading edge of the transition region (260 < X < 300) but not in 
the downstream or the so-called overshoot region where the magnetic field strength is the 
largest. The fact indicates that the observed energetic electrons are the result of reflection 
(not the leakage of downstream particles). The reflected electrons are accelerated parallel to 
the magnetic field during the reflection. 

In order to understand the electron acceleration process in more detail, individual tra- 
jectories of energetic electrons are analyzed. Figure H] shows a trajectory of typical energetic 
electron. The left panel represents the trajectory of electron (thick line) with the stacked 
profiles of the magnetic field B z (thin lines). The middle panel shows the time history of 
the particle energy. The solid, dotted and dashed lines of the panel display the perpendic- 
ular, parallel and total energy of the particle, respectively. The right-hand panel shows the 
time history of the first adiabatic invariant normalized to the upstream value. The unit of 
the vertical axis is the reciprocal of the electron gyrofrequency fi" 1 . The time history is 
plotted after Q ce t = 300 (u pe t = 6000). As is evident from the stack plot of the magnetic 
field profiles, the shock front is highly nonstationary and periodically reforms itself with the 
characteristic time scale of ~ 2fi^ i 1 . The shock wave is propagating toward the negative x 
direction with the average speed of ~ 0.5U . Initially, the particle located in the upstream 
is convected toward the shock with the E x B drift velocity. The particle first encounters 
the shock front at fl ce t ~ 70, and gains quickly its energy within a short time of VL ce At ~ 5. 
The energy of the particle increases by 2 orders of magnitude during the interval and the 
energy gain is almost perpendicular to the magnetic field. It is evident from the increase 
of the first adiabatic invariant that the process is highly nonad iabatic. This is indeed SSA 



mechanism investigated in detail by the previous studies (e.g. iHoshino fc Shimadal 120021 ). 
The properties of the acceleration are consistent with the previous results and we will not 
discuss the details of the process in any detail. After the first energization, the particle 
slowly drifts around the transition region without changing its energy until it collides with 
the overshoot at Q ce t ~ 180. The particle is pushed toward the upstream direction by the 
magnetic mirror force during the collision and gains its total energy, whereas the first adia- 
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batic invariant remains almost constant. In other words, the particle is reflected by the shock 
acting as a fast moving magnetic mirror and gains its momentum parallel to the magnetic 
field. The parallel energy increases, while the perpendicular energy is almost constant or 
slightly decreases during the interaction. 

The latter acceleration process associated with the reflection is known as SDA. As the 
shock angle and the upstream bulk flow speed increase, the energy gain of SDA increases. On 
the other hand, since the initial required energy for the reflection increases, the number of 
reflected particle rapidly decreases. Therefore, we expect that SDA process would not operate 
at high Mach number shocks. It is readily shown that the upstream thermal electrons of the 
present simulation cannot be reflected by the shock: The bulk velocities of the plasma in the 
upstream and downstream measured in the shock rest frame are U\/Uq ~ 1.5 and U2/U0 ~ 
0.5, respectively. The effective velocity of "the magnetic mirror" toward the upstream is 
U s /Uq = Ui/Uq/ cosObu — 8.6. The loss cone angle denoted by 9 C is given by 



where B\ and B max are the magnetic field strength in the upstream and the overshoot, 
respectively. The reflection will take place when the condition u± > U s sin 9 C is satisfied, 
where u± denotes the perpendicular velocity of particle measured in the upstream frame. 
If we take a typical magnetic compression ratio of B max /Bi ~ 10, the required velocity 
becomes U s sm6 c ~ 2.7Uq. This condition is quite severe, because the upstream electron 
thermal velocity is only 0.2f/ - Here, for simplicity, we ignore the effect of the electrostatic 
shock potential which suppresses the reflection efficiency. However, it cannot be neglected 
in general; hence, the reflection of thermal electrons requires even more stringent condition. 

In contrast to the above theoretical analysis without any nonadiabatic process, our 
simulation results demonstrate the reflection does indeed take place, which we attribute 
to the presence of the preacceleration via SSA. Figure [5] shows a schematic illustration of 
the reflection process initiated by SSA process. The cold upstream electrons are energized 
by the electrostatic waves excited by the Buneman instability at the leading edge of the 
shock transition region. The energization is so efficient that a nonnegligible fraction of 
electrons escapes outside the loss cone on a time scale of ~ O" 1 . The preaccelerated electrons 
escaping from the loss cone are subjected to SDA. In fact, by back tracing the trajectories 
of energetic electrons observed in the upstream side, we can confirm the scenario, i.e., the 
reflected electrons suffer the rapid energization via SSA when they enter the shock transition 
region and are reflected by the overshoot. On the other hand, low energy electrons are 
never reflected and just transmitted to the downstream. Therefore, we regard the two- 
step acceleration process as a preferential reflection process of energetic electrons. In the 
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Fig. 4. — Time history of typical energetic particle. From left to right, particle trajectory 
(thick line) and staked profiles of the magnetic field B z (thin lines), particle energy, the first 
adiabatic invariant, respectively. The solid, dotted and dashed lines of the middle panel 
show the perpendicular, parallel and total particle energy, respectively. Energy and the first 
adiabatic invariant are normalized to the upstream values. 
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present mechanism, SSA plays a crucial role to trigger the secondary acceleration process. 
The important point is that SSA produces suprathermal particles with approximately the 
power law energy spectra. Since the shock potential considerably reduces the reflection 
efficiency, the expected heating due to the Buneman instability is not sufficient to provide 
the required reflection energy and the production of suprathermal particles is essential. It is 
also important that the acceleration is almost perpendicular to the magnetic field. Because 
of this, pitch angles of energetic electrons become large and they are easily reflected by the 
shock. 



3. ELECTRON INJECTION MODEL 



The reflection of energetic electrons discussed in the previous section can be considered 
as "electron injection" to subsequent DSA process. The reflected electrons observed as a 
beam component in the upstream drive several instabilities. It is well known that fast electron 
stream parallel to the magnetic field excites Langmuir waves via the bump-on-tail instability. 
We observe the enhanced Langmuir turbulence in the upstream of the simulation results, 
which is not seen in strictly perpendicular shocks. Energetic electrons streaming away from 
the shock and associate d electrostatic turbulence are observed in the foreshock re gion of the 
Earth's bow shock (e.g. lAndersonlll969l ; I Anderson et al.lll98ll ; iKasaba et al.ll2000l ). Another 
electromagnetic instability may also be excited by electron cyclotron resonance. We expect a 
left-hand polarized electromagnetic wave propagating parallel to the beam which can scatter 
the energetic electron themselves, provided that the beam speed is greater than ~ va^-- As 
we will see in the discussion, this condition can be satisfied at shocks with Mach number 
typical of SNR shocks. Because of this, the initial field-aligned beam will tend to relax and 
become isotropic. From the above consideration, we recognize the reflected electrons as the 
seed population of DSA process, although a nonlinear evolution of the energetic electrons 
cannot be followed by our simulations because of the use of one-dimensional simulation box 
and a limited system size. 

The density of reflected electrons is identical to the injection efficiency defined as n^ T /n e , 
where n^ T is the number of nonthermal particles. Therefore, it is important to construct a 
theoretical model of the two-step acceleration process which predicts the injection efficiency 
and the energy density of nonthermal electrons relative to the thermal energy density. 
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Fig. 5. — Schematic illustration of surfing and drift acceleration. 
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3.1. Electron Heating and Acceleration in the Foot 



We investigate very high Mach number quasi-perpendicular shocks where the excitation 
of the Buneman instability is expected at the leading edge of the shock transition region. 
These shocks are typically observed around young SNRs. If we use the typical values of 
the magnetic field ~ 10/iG, the density ~ 0.1cm -3 and the shock speed ~ 10 4 km/s (e.g. 
Bamba et al.l 120031 ). the Mach number of the shock becomes Ma = 100 — 1000. Although 
temperatures of upstream plasmas are difficult to estimate, they are probably very cold 
{f3 e < 1). Therefore, we can expect that the threshold condition of the Buneman instability 



M A > JPe — 



(2) 



is satisfied at these shocks. 



In the following discussion, we neglect the relativistic effect for simplicity. However, 
we will show that the model agrees well with the simulation results which are in weakly 
relativistic regime (the maximum Lorentz factor of energetic electrons is 7 e ~ 2). 

Since the free energy of the Buneman instability is the relative streaming between the 
incoming electrons and the reflected ions, the amplitude of the electric field can be estimated 

as 



1 

:—i 
2 



a-n P m e V d ~ 



£2 
8vr' 



(3) 



where a is a energy conversion factor and Vd is the relative drift velocity (typically Vd ~ 2Vi). 
Although the precise value of a is hard to d etermine, it is a factor of the order of unity 
( llshihara et al.l Il98ll ; iDieckmann et al.l l2000bl ) . We assume that the wave electric field is 



parallel to the shock normal as in the case of our one-dimensional simulations. Thus, the 
electric fields parallel and perpendicular to the upstream magnetic field are respectively given 
by E\\ = E cos Obu and E± = E sin 9s n , where dsn is the shock angle. By equating the electric 
field energy to the electron thermal energy in each direction (parallel and perpendicular), 
we estimate the resultant electron thermal velocity v e ±, v e \\ after saturation as 



v e ,± — V\ sin 9 Bn 

V e \\ ~ Vi COS 6 Bn , 



(4) 
(5) 



respectively ( IPapadopoulosI 1 1 9881 ) . 



In addition to the strong electron heating, the simulation results clearly demonstrate the 
formation of high energy tail in the electron energy spectra due to the Buneman instability. 
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Therefore, we employ the bi-kappa distribution (jfJJ) as a model distribution function in the 
foot, 



f(v\\,v±) 



n .f 



r(« + i) 



,±^e,|| (7Tk) 3/2 T(k- 1/2) 



1 i (vj (n-Vs h ) 



K \ V~ 



-K-l 



(6) 



where is the gamma function, V s h and nj are the parallel drift velocity and the density 
in the foot, respectively. The distribution function is measured in the HTF. We define 
the foot as the region where the Buneman instability saturates. Hereafter, the subscript / 
represents the value in the foot. The incoming electrons are decelerated so as to cancel the 
zeroth order current, and somewhat compressed (rif > The mass conservation law leads 
to Vf = ViUi/nf. The magnetic field is also compressed with the same compression ratio 
Bf = Biiif/nx- Thus, the parallel drift velocity of the compressed incoming electron is given 
by 



v sh 



cos 9 



(7) 



Brn 



which is almost equal to V\j cos9 Bn for 9 Bn ~ 90°. Here 6b fn is defined as the angle between 
the shock normal and the magnetic field line in the foot. 



3.2. Adiabatic Reflection 

After the energization due to the Buneman instability, the ion acoustic instability is 
triggered and further electron heating may occur. However, this heating is not important 
for the generation of suprathermal electrons. Since we find that the first adiabatic invariant 
of energetic electron is conserved during the reflection process, we use the adiabatic approx- 
imation after the energization due to the Buneman instability. This approximation enables 
us to compute the reflected electron density by integrating the distribution function outside 
the loss cone. In general, the electrostatic shock potential is known to affect the mirror 
reflection process. If we include the finite potential measured in the HTF <f) HT , the condition 
of the electron reflection becomes 

v l> (v?\+—<f> Err )taiie e (8) 

V m e J 

where 6 C is the loss cone angle given by 



(9) 
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Here we use the foot to the maximum compression ratio to estimate the loss cone angle, 
because we use the adiabatic approximation after the energization in the foot. We consider 
Bmax as the magnetic field strength in the overshoot. We know from both observations and 
numerical simulations that the shock potential is determined by the upstream bulk ion flow 
energy. Therefore, we normalize the shock potential to the bulk ion energy as 



lHT 



hrnVl 



2 ' 



(10) 



The effect of the potential cannot be neglected when the particle parallel energy is comparable 
or smaller than the potential. Since the typical parallel velocity is expressed as v\\ ~ V s h — 
V\j cos 9 Bni the condition can be written as 



<p 



COS 2 9 Bn > 1. 



For (f) = 0.4 and rrii/m e = 100, this condition leads to 9 Bn < 81°. Moreover, this critical 
angle increases with increasing mijm e . The use of the realistic proton to electron mass ratio 
gives 9bu 88°. Consequently, the effect of the shock potential is important for a wide 
range of shock angles and should be included in the injection model. We should note that 
the use of unrealistically small value of mj/m e , which is usual in most of PIC simulations, 
overestimates the reflection efficiency. 

Although the shock potential in the HTF is identical to that in the NIF in the absence 
of the noncoplanar magnetic field component, significant field rotations out of the copla- 
narity plane within the shock tr ansition region are of ten found by both in situ observations 



and numerical simulations (e.g. iThomsen et all 119871 ). Our numerical simulations also pro- 



duce the noncoplanar magn etic field component (B y ) in the transition region. However, 
Krauss-Varban &: Wul (119891 ) showed that the effect of the noncoplanar magnetic field com- 
ponent is only of second order in B y /B and thus has little influence on the electron kinetics. 
Therefore, we can use the potential measured in the NIF as <p HT . 



The reflected electron density n r can be written as 



n r 



2tt / dv\\ 
'o 



„2 _l- 2^^HT 



v±dvxf(v\\,v ± ) 



(12) 



tan t 



Evaluating the integral analytically, we obtain 



n 



11 ,. 



f 



2v P \\r K \/KTi 




2q T(k) 



pT(/c-l/2j l 2' K '2' rp J 



(13) 
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where F(a,(3, 7; x) is the Gauss hypergeometric function and p, q, and r are respectively 
given by 

r + tan 2 # c 

p = (14) 

9 = ^ (15) 

■^-0 HT tan 2 c vl rtan 2 # 

r = ^-— 2 -+ V f ^ 2 C Q + 1, (16) 

/tf g j_ /tf g j_ r + tan 2 C 

using the anisotropy defined as r = v\^l v \ 11 ■ 

The reflected electron beam velocity measured in the shock frame can be easily estimated 
from the adiabatic theory as 



T r \ B max Bf 
Vr = Vsh \j B f - Bl ■ (17) 

The beam temperature can be determined by the integral 

/* 00 n 00 J 

n r T r = 2ir dv\\ v±dv±-m e (v» + vj_)f(v\\, v±). (18) 

Jo y^„|+^^HT tanflc 2 

However, we use a rather simple estimate instead of this. The thermal energy density of the 
reflected part of the distribution function can be written as 

i>co poo 

E th = / N(E)EdE = / AE~ k dE, (19) 
Je Je 

where A and E min are the normalization constant and the minimum energy required for the 
reflection, respectively. In the above integral, we approximate the energy spectrum by the 
power law part of the distribution function and the anisotropy is neglected. Substituting the 
normalization constant A = n r KE^ in , we obtain 

E th = n r ^-E miQ . (20) 

K — 1 

The minimum energy is approximated by 

£ min ~ Qm e V^ + e<p HT ^j tan 2 9 C . (21) 

By combining these equations, we can write the energy density of the beam e r as 

1 k f\ \ 

e r ~ -n r m e vl H -n r -m e Vf h + e<f) HT tan 2 9 C . (22) 

2 K — 1 \Z J 
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This is somewhat crude estimate. Since the reflected beam distribution will be modified 
by nonadiabatic processes in the shock transition region, even the energy density estimated 
from the analytic integration (ITS]) might be inaccurate. Nonetheless, (122]) gives reasonably 
good estimate if it is compared to the simulation results. Therefore, we adopt this simple 
approximation. In the following discussion, density, velocity and energy will be given in 
units of the density, bulk velocity and bulk electron energy in the upstream, respectively. 
Note that, under this normalization, the result of the present model is independent of M4 
and @ e , provided that ([2J is satisfied. 



3.3. Comparison Between Model and Simulation Results 

Obviously, the injection efficiency strongly depends on the choice of k. We adopt k = 2.5 
(corresponding to the power law index of 3.5) which approximates the energy spectra in the 
shock transition region (see Figure [SJ. Likewise, the value of = 0.4 is used for the model 
calculation. The results depend weakly on the magnetic field in the foot and the overshoot. 
In the following discussion, we use the fixed values of Bf/B\ = 2.5 and B max /Bi = 10.0 
which are the typical values obtained from the simulation results. The variation of these 
values does not affect the result significantly. 

Figure O shows the shock angle dependence of the density and the energy density of the 
reflected electrons calculated by the model with those obtained from the simulation results. 
The energy density of the model is converted to the downstream rest frame. The simulation 
results are averaged over X s h — 70c/u pe < x < X s h — 20c/ 'u pe , where X s h is the position 
at which the magnetic field is compressed by a factor of 1.5 from the upstream value. The 
time-averaged values during the last interval AT ~ fi" 1 of each run are shown. The error 
bars correspond to the variance during the interval. 

The model curves agree well with the simulation results when the shock angle is small 
(#Bn ^ 75). However, the simulation results tend to deviate from the model curves with 
increasing the shock angle. At the shock angle of Qbu — 85°, no reflected electrons are 
observed in the simulation, while the model gives the density of ~ 4 x 10~ 3 . This large 
difference is not acceptable when we apply the model to real shocks. In the next section, we 
introduce some modifications into the model in order to correct the discrepancy. 
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Fig. 6. — Comparison between our theoretical model and the simulation results. The solid 
and the dashed lines show density and energy density of reflected electrons calculated by the 
model. The simulation results are displayed by filled circles (density) and crosses (energy 
density) . 
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3.4. Corrective Effects 

3.4-1- Maximum Energy Cut-Off of Shock Surfing Acceleration 

The actual energy spectra obtained from the simulations have a cut-off energy where 
the number of particle rapidly falls off. Therefore, the number of reflected electron will fall 
off as oc exp(—E min /E max ), where E max is the maximum cut-off energy of SSA. It is easy 
to understand that the presence of the maximum cut-off energy becomes important when 
it is comparable or smaller than E min . Because E min increases with the shock angle, the 
presence of the maximum energy introduces the cut-off shock angle above which the reflection 
efficiency rapidly falls off. We can determine E max ~ 100 from the simulation results. This 
leads to the cut-off shock angle of ~ 86°, which is larger than that of the simulation results. 

Although only the above correction effect cannot explain why the reflected electron 
density of the simulation results show the rapid decrease at the shock angle #g n > 80, the 
effect may be still important even when the shock angle is smaller than the cut-off, if we 
consider the application to real shocks. This is due to the fact that E min strongly depends 
on the shock potential which increases with increasing the mass ratio. As discussed later, 
E max should be larger than the shock potential in order to obtain a large fraction of reflected 
electrons. 



3-4-2. Escape Probability During Reflection 

Although we use the adiabatic approximation after the first energization to compute 
the reflected electron density, the first adiabatic invariant may be violated in the presence of 
nonstationarity of the shock structure and/or the ion acoustic turbulence in the transition 
region. Particles which are initially outside the loss cone may be scattered by these effects 
and fall inside the loss cone. If we consider the escape probability P esc as a constant both in 
time and space during the reflection, the number density will fall off as oc exp (— P esc T re f), 
where T re f is the characteristic time for the reflection process. We can write T re f by using 
initial perpendicular particle energy E re f (see Appendix) 

o T , 2m e V sl (2 o\ 

"o'ref- Eref / tfm 2 dc _ e(f> HT> V 6 > 

where the spatial gradient of the magnetic field strength and the electrostatic potential are 
assumed to be constant and the shock width is evaluated as Vi/Q C i- We approximate the 
typical reflected electron energy before the reflection process as E re f ~ (E max + E min )/2. 
It is readily shown that T re f increases with increasing the shock angle, which leads to the 
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reduction of the reflection efficiency. We should note that the parameter P esc is introduced 
as a free parameter to correct the discrepancy between the model and the simulation results. 
We determine the value of P esc so that the model fits the simulation results. 

Thereafter, P esc will be given in the units of ion gyrofrequency Q C i, because T re f is of 
the order of ion gyroperiod. 



Rewriting the density given by (IT51) as n' r , above two effects can be included in the 
model with the following form 



Figure [7] shows a comparison between the simulation results and the model using the 
above correction. The format of the figure is the same as that of Figure [6j We use P esc = 4 
and E max = 100 to produce the figure. In this case, we see that the model is significantly 
improved and the agreement becomes quite good. Therefore, we can conclude that the 
present model can predict the injection efficiency as well as the energy density of nonthermal 
particles, if the model parameters are properly estimated. 

Note that the value P esc = 4 used here indicates that the underlying mechanism of elec- 
tron scattering in the shock transition region is the self- reformation of the shock front, which 
occurs with the characteristic time scale of ~ 2S7 ( I^ 1 . Since the shock structure considerably 
changes with this time scale, the motion of electrons interacting with the shock longer than 
this period may also be modified. As a result, only electrons which have sufficient energy 
for the rapid reflection can escape to the upstream. 



Now, let us apply the model to SNR shocks. First, we should discuss the maximum 
energy of SSA. Since the shock potential is determined by the ion bulk flow energy, it will 
increase with increasing the mass ratio. The use of the same maximum energy E max for the 
real mass ratio shock will result in the net reduction of the injection efficiency. 



3.4-3. Comparison Between Corrected Model and Simulation Results 
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3.5. Application to Supernova Shocks 
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Fig. 7. — Comparison between our theoretical model with corrections and the simulation 
results. The format is the same as Figure [6j P esc = 4 and E max = 100 are used to produce 
the figure. 
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Hoshino &: Shimadal (120021 ) discussed the maximum energy of SSA and obtained the 



condition 



/ 771 

2M A Ja—>\ (25) 

V m 

for "unlimited electron acceleration", which means trapped electrons cannot escape from 
ESWs and continue to accelerate. They also mentioned that even in this regime, the trapping 
time might be limited by some other important factors. These include shock front nonstation- 
arity (e.g. self-reformation) and multidimensional effects. In addit ion to this, the stability of 



ESW s is also a quite important factor limiting the trapping time. iDieckmann et al.l ((2000b 



20041 ) investigated the lifetime of the BGK mode excited by the Buneman instability using 
a periodic simulation model. They showed that the lifetime depends on the drift velocity of 
the reflected ions and also the background ion temperature. In our curr ent understanding 



i t is s till a hard task to estimate the trapping time scale. Furthermore, iMcClements et al. 



(120011 ) showed that stochastic electron acceleration can occur when u pe /VL ce 3> 1. In other 
words, electrons once detrapped from the potential can interact with the wave again. In 
this case, the maximum energy of SSA is not simply limited by the trapping time scale. It 
is very difficult to estimate the maximum energy of SSA and far beyond the scope of the 
present paper. We simply assume that the maximum energy increases with increasing the 
mass ratio from the following qualitative consideration: If the mass ratio is increased, the 
free energy provided by the reflected ions also increases. Since the free energy is the source 
of the energetic electrons, the maximum energy of SSA will also become higher. 

If the maximum energy is smaller than the shock potential, the injection efficiency is 
greatly reduced independent of the shock angle. For typical values of the shock potential 
~ 0.4, we require E max > 750 in order to obtain a measurable fraction of the reflected 
electrons. We assume E max = 1000 in the following discussion, which leads to the cut-off 
shock angle of ~ 88°. In this case, the presence of the maximum energy cut-off does not 
significantly affect the result when the shock angle is smaller than the cut-off shock angle. 
Note that E max = 1000 corresponds to 10 - 100 keV at typical SNR shocks (Vi/c ~ 10~ 2 ). 

The escape probability P esc is also important in the sense that it introduces another 
cut-off shock angle. If P esc is determined by the frequency of the self-reformation, it should 
be independent of the mass ratio. Therefore, we adopt P esc = 4 obtained by the simulation 
results of mj/m e = 100. This leads to the cut-off shock angle of ~ 85°, which is smaller than 
that introduced by E max . Namely, in the present parameter range, P esc is more important 
than E max for determining the shock angle dependence. 

Figure [8] shows the injection efficiency (left) and the energy density (right) obtained by 
the present model. We consider the shock potential as a free parameter, because it is also 
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difficult to estimate the precise value of the shock potential at high Mach number shocks. 
We expect the potential will be in the range of (f) = 0.3 — 0.5. Then, we obtain the peak 
injection efficiency of ~ 2 x 10~ 4 at Obu — 80°. Similarly, the energy density at the peak 
shock angle is approximately 10%. 

These re s ults can be directly compared to X-ray and radio observations of SNRs. 
Bamba et al.l (120031 ) carried out a detailed investigation of the northeast shell of SN 1006 
observed by Chandra and argued that the estimated injection efficiency is ~ 1 x 10~ 3 and the 
energy density of nonthermal particles is about 30% of the therm al energy density. These 
values are similar to previous observations (e.g. lAllen et al.ll200ll ). The model prediction 



shows a good agreement with the observation, although the observed injection efficiency 
and the energy density of cosmic ray electrons are slightly larger than those obtained by 
the present model. Furthermore, the high resolution observation reveals that nonthermal 
emission is confined in very thin filaments. The spatial scales of the nonthermal filaments 
are ~ 0.04pc an d ~ 0.2pc in the up stream and downstream, respectively. Similar results are 
also reported by lLong et al.l (120031 ) . If the standard DSA theory is employed, the observed 
scale length of the nonthermal filaments imposes the constraint that the shock angle should 
be sufficiently large [Qbu ^ 80), unless we assume the m agnetic field of 20 — 85/xG, which is 



large r than the usual interstellar value of a few /zG (e.g. iBamba et al.ll2003l ; lYamazaki et al. 
2004 ). The peak shock angle predicted by the model is around 80°, again we see a good 
agreement with the observation. 

We should note that the peak shock angle depends on the choice of both P esc and 
Emax- If E max is sufficiently large compared to the shock potential, the peak shock angle is 
simply determined by P esc . Although we attribute the escape probability to the shock front 
nonstationarity in the present analysis, there remain some other possibilities such as pitch 
angle scattering due to the interaction with the whistler waves in the shock transition region. 
The discussion of such nonadiabatic behaviors of energetic particles requires more detailed 
understandings of the structure and the wave activities in the shock transition region. 



4. DISCUSSION 

We have studied rapid electron energization mechanism within the transition region of 
high Mach number quasi-perpendicular shocks. We found that highly energetic electrons are 
generated through successive two different acceleration processes. First, energetic electrons 
are produced via SSA at the leading edge of the shock transition region. As a result, the 
preaccelerated electrons escaping outside the loss cone are subjected to SDA and preferen- 
tially reflected back to the upstream. We consider the two-step acceleration mechanism as 
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the injection to subsequent DSA process. We have constructed a model of the acceleration 
mechanism which predicts the injection efficiency and the energy density of nonthermal par- 
ticles. The estimated injection efficiency agrees well with observations of SN 1006. We also 
found that the shock angle dependence of the injection efficiency is consistent with the shock 
angle constraint inferred from observations. 

Although the present model generally agrees well with observations, there remain some 
important issues. The most important one is the acceleration efficiency of SSA, i.e. the 
maximum energy and the spectral index. It is easily understood that the spectral index 
affect the injection efficiency to a great extent. We have used the power law index of 3.5 
throughout in this paper. However, it may depend on some important physical parameters 
such as uj pe /Q ce , #e n , mi/m e etc. We also observe that the spectral index varies with time 
according to the phase of the self-reformation. The maximum energy of SSA is also important 
for the electron injection. The maximum energy, which is required in order to account for 
observations, depends on the shock potential. The maximum energy should be larger than 
the shock potential to obtain a measurable fraction of reflected electrons. Even in this 
case, the maximum energy is important for determining the shock angle dependence of the 
injection efficiency, because it introduces the cut-off shock angle where the injection efficiency 
rapidly decreases. 

We also introduce the escape probability of energetic electrons as another important 
factor for determining the shock angle dependence. We attribute the probability to the self- 
reformation of the shock front. This is because the value of the escape probability estimated 
from the simulation results is P esc ~ 4fi C j, which indicates that the escape mechanism is, 
to some extent, related to the shock self- reformation process. Although the self -reformation 



process of the shock fron t has been extensive ly studied by many authors (e.g. lQuestlll986 



Lembege fc Savoinil 1992 ; Scholer et a 



2003), it is still a controversial topic of collisionless 



shock physics. IShimada fc Hoshinol (120051 ) have recently pointed out that the dynamics 



of the shock front can be modified by strong dissipation due to microinstabilities in the 
transition region of very high Mach numbe r shocks. Similar phenomen on at moderate Mach 
number shock Ma ~ 6 is also reported by IScholer fe Matsukivol (120041). wher e the modified 
two stream instability plays an important role (e.g. lMatsukiyo fe Scholerll2003l ). It is difficult 
to say whether the shock self-reformation process survives in real SNR shocks, where strong 
dissipation of both electrons and ions is expected. If the self-reformation is suppressed by 
strong dissipation, the escape probability may become smaller. If this is the case, the cut- 
off shock angle determined by the escape probability disappears and the maximum cut-off 
energy of SSA controls the shock angle dependence. 



If the shock self-reformation is suppressed and the maximum energy of SSA becomes 
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sufficiently large, the reflection of energetic electrons will take place at the shock angle 
very close to the threshold between subluminal and superluminal shocks. The relativistic 
effect becomes important for such situations, because both the effective shock velocity V s h 
and the required energy for the reflection become relativistic. In this case, since the shock 
angle is well close to 90°, the shock potential has little influence on the reflection process, 
i.e., the relativistic effect does not affect the mirror reflection efficiency in the absence of 
the electrostatic potential. Thus, the present mechanism should work in principle. It is 
important to investigate whether highly relativistic electrons are generat ed via SSA or not. 



Stoch astic version of SSA may be important to understand the issue (IMcClements et al. 



200lh . 



We should also point out that further electron heating and acceleration in the shock 
transition region may be possible if multidimensional effects are considered. It is well known 
that many plas ma microinstabilities can be excited within the shock transition region (e.g. 



Wu et al.lll984l ). For instance, it is easy to expect that the whistler mode waves are excited 
by temperature anisotropy (T± > T\\) of thermal electron because of strong perpendicular 
heating and acceleration due to SSA and the ion acoustic turbulence. We also expect that the 
reflected electron beam with temperature anisotropy (T± > Tji) will excite the w histler waves 



propa gating antiparallel to the beam by cyclotron resonance as discussed by iTokar et al. 



( 119841 ). although the use of one-dimensionality assumption in the present simulation inhibits 
the excitation of these instabilities. If the whistler wave intensity becomes sufficiently strong, 
both the beam and core electron distributions will become isotropic via strong pitch angle 
scattering. The first adiabatic invariant of electron will be violated and our simple theoretical 
model using the adiabatic approximation may be inaccurate. Even in this case, there are 
no reason why SDA should not operate, because the physical mechanism is quite simple and 
does not require any special conditions. However, the process will be strongly modified by 
the turbulence in the shock transition region. The wave-particle interaction within the shock 
layer may provide further heating and acceleration of electrons. We would like to emphasize 
again the significant importance of SSA on the energization of electrons. SSA plays a key 
role in the turbulent shock structure in the sense that it provides additional sources of free 
energy and may lead to further energization of electrons. 

The relation between th e elect ron acceleration efficiency and the whistler waves is re- 



cently studied by lOka et al.l (120061 ). They analyzed a number of the Earth's bow shock 
crossing events observed by Geotail. They clearly showed that the power law index of elec- 
tron energy spectra measured in the shock transition region is regulated by the so-called 
whistler critical Mach number M™ it , which is defined as the critical point above which the 
whistler waves cannot propagate upstream. The spectral indices are distributed 3.5 — 5.0 in 
the sub-critical regime, while the harder energy spectra with indices of 3.0 — 3.5 are observed 
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in the super-critical regime. In the super-critical regime, the whistler waves generated by 
microinstabilities are accumulated in the shock transition region. Since SNR shocks are in 
super-critical regime, we can expect that the accumulated energy of the whistler waves may 
contribute to further electron energization. It is interesting to investigate the relationship 
between the whistler wa ves and the electron accelerat i on ef ficiency. Gyroresonant surfing 



acceleration proposed by lKuramitsu fc Krasnoselskikhl (120051 ) might be important. 



In order to model such nonadiabatic processes, we must know the shock structure (e.g. 
shock potential, nonstationarity) and the wave activity in the transition region of realis- 
tic high Mach number shocks in more detail. Numerical simulation of self-consistent shock 
structures including the whistler wave turbulence requires at least two-dimensional simu- 
lation domain which demands very large computational resources. Another possibility to 
improve the understandings of high Mach number shocks is in situ observation of interplane- 
tary shocks in the inner heliosphere. It is known that interplanetary shocks dri ven by Coronal 



Mass Ejections (CMEs) can be very high Mach number near the Sun (e.g. ISmart fc Shea 



19851 ). Expected Mach number at the Mercury orbit of ~ 0.4 AU becomes Ma ~ 40. Obser- 
vations of such interplanetary shocks may provide us useful information to understand the 
physics of electron acceleration at very high Mach number shocks. 

In the present paper, we have restricted ourselves to the discussion of the electron 
injection process. The self-generation of upstream waves is one of the major problems of 
the electron DSA theory. The cyclotron resonance condition between the reflected electron 
beam and a left-hand polarized weakly damped Alfven wave (/ci^/f2 C j < 1) requires, 

Vr- > l_J_nH (26) 

where the factor 1/2 in the right-hand side of the above equation indicates the difference of 
the frame between the shock frame and the upstream frame. If we use typical Mach number 
of SNR shocks Ma ~ 100 — 1000, this condition becomes v r /V\ > 1 — 10. By combining 
ffT7|) and (1261) . we can estimate the critical Mach number above which the self-generation of 
upstream waves becomes possible. Figure [9] shows the shock angle dependence of the critical 
Mach number. We confirm that the self-generation of upstream waves is indeed possible in 
typical SNR shocks, i.e., our theory successfully explains the electron injection as well as the 
triggering of subsequent DSA process. 

Since we also know the reflected electron density, we can, in principle, estimate the wave 
intensity, hence, the diffusion coefficient of energetic electrons using quasi-linear theory. 
However, we think that it is rather important to investigate further electron energization 
in the shock transition region introduced by multidimensional effect before we discuss the 
consequence of the electron injection and subsequent DSA process. 
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Another important problem is the back-reaction from nonthermal particles. It is known 
that the energetic particles affect the upstream plasma environme nt when their energy den- 



sity becomes comparable to that of the background plasma (e.g. iDrury fc Volk!ll98ll ). We 
cannot discuss such a nonlinear evolution in the present model, because we do not consider 
the injection process of ions. Obviously, ions have much larger energy density than that of 
electrons. Thus, the injection efficiency of ions is more important to understand the shock 
structure in the presence of energetic particles. In order to study the nonlinear evolution in- 
cluding the interaction between thermal and nonthermal particles, the injection efficiencies 
of both electrons and ions will be of great importance. Understanding of both the injec- 
tion processes and the nonlinear shock structure will elucidate the problem of cosmic ray 
acceleration at collisionless shock waves. 
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A. Calculation of Reflection Time 

The parallel equation of motion in the HTF under the action of the magnetic mirror 
force and the electrostatic potential can be written as 

dv\\ ( dB d<f>\ 

m e-J7 = ~ + ^ COS 9 Bn , (Al) 



dt \ dx dx r 

where x represents the particle position and \x = m e v\/2B is the first adiabatic invariant. 
By assuming the spatial gradients of both the magnetic field strength and the electrostatic 
potential are constant, we can integrate the equation, 

w 11 w 2m e Ax K J 

In this calculation, we have neglected the curvature of the magnetic field line in the shock 
transition region for simplicity (i.e. cos 9 Bn = const.). The initial parallel velocity is v\\(0) ~ 
V s h and the shock width is approximately given by Ax ~ Vi/Vt ci . Then, we obtain the 
characteristic time for the reflection process T re f as 



(A3) 
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If we use the definition of the characteristic particle energy perpendicular to the magnetic 
field before the reflection E re f = m e Vj_/2, 



/iAi? = E re f 



B max — Bf 



E re f 

tan 2 6 C 

Thus, we can rewrite the reflection time 



Bf 

(A4) 



2m V 2 

Q aT ~ , f ah — . (A5) 

f £ re/ /tan 2 # c -eAc/> V ; 

It is easy to understand that T re f is proportional to 1/ cos 2 Osn f° r the fixed particle energy. 
The reflection time becomes shorter with increasing the particle energy. 
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Fig. 8. — Density (left) and energy density (right) of reflected electrons. The contour levels 
are equally spaced on a logarithmic scale. The vertical axis represents the shock potential 
normalized to the ion bulk energy in the upstream. 
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Fig. 9. — Critical Mach number above which the self-generation of upstream waves becomes 
possible. 



